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1 Introduction 



The apparent success of the standard model (SM) in describing the elementary particles 
and their interactions is spoiled by the presence of the plethora of free parameters. There- 
fore, it provides a challenge to theorists to attempt to understand at least some of these 
parameters. 

The traditional way to reduce the independent parameters of a theory is the introduc- 
tion of a symmetry. Grand Unified Theories (GUTs) |], 0, |3|] are representative examples 
of such attempts, and some of them are certainly successful to some extent. For instance, 
the minimal SU (5) reduces by one the gauge couplings of the SM and gives us a testable 
prediction for one of them. In fact, LEP data seem to suggest that a further symmetry, 
namely N = 1 global supersymmetry [|5| should also be required to make the prediction 
viable. Furthermore, GUTs also relate Yukawa couplings among themselves, which in 
turn might lead to testable predictions for the parameters of the SM. The prediction of 
the ratio m T /mb |6| in the minimal SU(5) was an example of a successful reduction of the 
independent parameters of this sector of the theory. On the other hand, requiring more 
symmetry (e.g 5*0(10), E 6 , E 7 , E 8 ) does not necessarily lead to more predictions for the 
SM parameters, due to the presence of new degrees of freedom, various ways and channels 
of breaking the theory, etc. An extreme case from this point of view are superstrings, 
which have huge symmetries, but no predictions for the SM parameters. 

In a series of papers J7], [5[ |TD[ we have proposed that a natural gradual extension 
of the GUTs ideas, which preserves their successes and enhances the predictions, is to 
attempt to relate the gauge and Yukawa couplings of a GUT, or in other words, to 
achieve Gauge- Yukawa Unification (GYU). 

Searching for a symmetry that could provide such a unification, one is led to introduce 
a symmetry that relates fields with different spins, i.e. supersymmetry, and in particular 
N = 2 supersymmetry |TJ]. Unfortunately N = 2 supersymmetric theories have serious 
phenomenological problems due to light mirror fermions. We expect that a GYU is a 
functional relationship which is derived by some principle. In superstring theories or in 
composite models there exist relations among gauge and Yukawa couplings, although in 
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practice both kinds of theories have more problems than the SM. 

Before turning to our attempts, let us mention some earlier ones in order to make 



clear which are really the predictions in each case. Veltman, Decker and Pestieau |L2 
by requiring the absence of quadratic divergences in the SM, found that the following 
relationship has to hold: 

m e + m fi + m r + ^( m u + m d + m c + m s + m t + m b) 

3 o 3 o 3 9 /-i \ 

= + l m z + 1 m H I 1 ) 

A very similar relation is obtained demanding spontaneous breaking of supersymmetry 
via F-terms. In that case one obtains 1H 



E(-1) 2J (2J + 1H 2 = 0, (2) 

3 

where j is the spin of the particle. When this formula is applied to the SM we obtain a 
relation which differs from Eq.(|l|) only in the coefficient of m? H . In both cases a prediction 
for the top quark was possible only when it was permitted experimentally to neglect 
the ran as compared to mw,z with the result m t = 69 GeV. Otherwise there is only a 
quadratic relation among m t and m#. 

A celebrated relation among gauge and Yukawa couplings is the Pendleton-Ross (P- 
R) infrared fixed point [|14| . The P-R proposal, involving the Yukawa coupling of the top 
quark g t and the strong gauge coupling as, was that the ratio a t /as, where a t = g^f^n, 
has an infrared fixed point. This assumption predicted m t ~ 100 GeV, and therefore it 
is ruled out. Moreover, it has been shown that the P-R conjecture is not justified at 
two- loops. On the contrary, the ratio at/ diverges in the infrared. 

Another interesting conjecture, made by Hill [I5|, is that at itself develops an infrared 



fixed point, leading to the prediction m t ~ 280 GeV. 

The P-R and Hill conjectures have been done in the framework on the SM. The 
same conjectures within the minimal supersymmetric SM (MSSM) lead to the following 
relations: 

m t ~ 140 GeV sin (3 (P-R) (3) 
m t ~ 200 GeV sin (3 (Hill) (4) 



where tan/3 = v u /vd is the ratio of the two VEV of the Higgs fields of the MSSM. We 
should stress that in this case there is no prediction for m t , given that sin/3 is not fixed 
from other considerations. Therefore, the conclusion is that all the attempts that have 
been made so far to extract predictions by some principle, leading to relations among 
gauge and Yukawa couplings were found either wrong or not predictive enough. 

In the following we would like to emphasize an alternative way to achieve unification of 
couplings, which is based on the fact that within the framework of a renormalizable field 
theory, one can find renormalization group invariant (RGI) relations among parameters, 
that can improve the calculability and the predictive power of a theory. In our recent 
studies [0, ||, |9|, [HJ, we have considered the GYU which is based on the principles of 
reduction of couplings |TB|, [T7|, JSJ, |TS], fZIJ and finiteness [g f^, 0, f?§ [25|, |2§, |27 



These principles, which are formulated in perturbation theory, are not explicit symmetry 
principles, although they might imply symmetries. The former principle is based on the 
existence of RGI relations among couplings, which preserve perturbative renormalizability. 
Similarly, the latter one is based on the fact that it is possible to find RGI relations among 
couplings that keep finiteness in perturbation theory, even to all orders. Applying these 
principles one can relate the gauge and Yukawa couplings without introducing necessarily 
a symmetry, nevertheless improving the predictive power of a model. 

2 Unification of Couplings by the RGI Method 

Let us next briefly outline the idea. Any RGI relation among couplings (which does not 
depend on the renormalization scale fi explicitly) can be expressed, in the implicit form 
^(^l) ' " ' j 9a) = const., which has to satisfy the partial differential equation (PDE) 

d{i S og a 

where j3 a is the /3-function of g a . This PDE is equivalent to the set to ordinary differential 
equations, the so-called reduction equations (REs) flT7| , 



A? -j 1 = @ a ' a = 1 > ' ' ' ' A ' ( 6 ) 

dg 
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where g and (3 g are the primary coupling and its /^-function, and a does not include it. 
Since maximally (A — 1) independent RGI "constraints" in the A-dimensional space of 
couplings can be imposed by $ a 's, one could in principle express all the couplings in terms 
of a single coupling g. The strongest requirement is to demand power series solutions to 
the REs, 

9a = Epi" + V n+1 , (T) 

n=0 

which formally preserve perturbative renormalizability. Remarkably, the uniqueness of 
such power series solutions can be decided already at the one-loop level |L7] . To illustrate 
this, let us assume that the /^-functions have the form 



p a = tt^ [E Pa 1)bcd g»gcg d + Y,A 1)h g b g 2 ] + 

b,c,dytg b^g 



where • • • stands for higher order terms, and (3^ bcd 's are symmetric in b,c,d. We then 
assume that the p^'s with n < r have been uniquely determined. To obtain p( r+1 )'s, we 
insert the power series into the REs (^) and collect terms of 0(g 2r+3 ) and find 

E M(r)a p d r = lower order quantities , 

d+g 

where the r.h.s. is known by assumption, and 

M{r) d a = 3 E bcd Pb ] P ( c ] + f3 a 1] d -(2r + l) /3« 5 d a , (9) 

o = E^^^fVM' + Etf^f'-^pi 1 '. (io) 

b,c,d^g d^g 

Therefore, the p^'s for all n > 1 for a given set of p^'s can be uniquely determined if 
det M(n) d a ^ for all n > 0. 

The possibility of coupling unification described above is without any doubt attractive 
because the "completely reduced" theory contains only one independent coupling, but it 
can be unrealistic. Therefore, one often would like to impose fewer RGI constraints, and 



this is the idea of partial reduction []TJ 



In the following chapters we would like to consider three different GYU models based 
on supersymmetric unified theories. Before doing so let us recall here some features 



which are common in all cases. Let us consider a chiral, anomaly free, N — 1 globally 
supersymmetric gauge theory based on a group G with the gauge coupling constant g. 
The superpotential of the theory is given by 

W = ^rriij^i^j + ^C ijk <j)i<f) j <f) k , (11) 

where rriij and Cij k are gauge invariant tensors and the matter field fa transforms according 
to the irreducible representation Ri of the gauge group G. The renormalization constants 
associated with the superpotential fllTD , assuming that supersymmetry is preserved, are 



0° = (^) (1/2) 0,, 


(12) 


ryi'j' 
ll'ij = Z ij 111 -'J' , 


(13) 


n Q yi'j'k' n 
^ijk — ^ijk ^i'j'k 1 ■ 


(14) 



The N = 1 non-renormalization theorem |30| ensures that there are no mass and cubic- 
interaction-term infinities and therefore 

yi'j'k' 7 l/2i" „l/2j» „l/2k" _ ri"ri"rfe" 
^ijk ^i' ^j' ~ °(i °j °k) > 

Zifzir'zp" = 8% 8%. (15) 

As a result the only surviving possible infinities are the wave-function renormalization 
constants Zf, i.e., one infinity for each field. The one -loop /^-function of the gauge 
coupling g is given by E3| 



* ) = l = i&E w- 3 ^]. ( l6 ) 

where l{Rj) is the Dynkin index of Ri and C^iG) is the quadratic Casimir of the adjoint 
representation of the gauge group G. The /^-functions of Cijk, by virtue of the non- 
renormalization theorem, are related to the anomalous dimension matrix 7^ of the matter 
fields fa as: 

frjk = = Ciji j l k + C lkl ^ + C jkl 7 j . (17) 



At one-loop level 7^ is |22 



where C^-Rj) is the quadratic Casimir of the representation Ri, and C l ^ k = C*j k . Since 
dimensional coupling parameters such as masses and couplings of cubic scalar field terms 
do not influence the asymptotic properties of a theory on which we are interested here, it is 
sufficient to take into account only the dimensionless supersymmetric couplings such as g 
and Cijk- So we neglect the existence of dimensional parameters, and assume furthermore 
that Cij]^ are real so that Cfj k always are positive numbers. For our purposes, it is 
convenient to work with the square of the couplings and to arrange Cij k in such a way 
that they are covered by a single index % (i = 1, • • • , n): 

a = 4^ ' ai = 17 ■ (19) 
The evolution equations of a's in perturbation theory then take the form 

Tt = P = -^ (1) « 2 + --- > 

= A = -ti 1) a i a + '£/3Q k a j a k + --- , (20) 

where • ■ ■ denotes the contributions from higher orders, and (3^ k = Ajy- 

Given the set of the evolution equations (12(1), we investigate the asymptotic properties, 



as follows. First we define 16, 17 



Si = — , i = 1, • • • ,n , (21) 



a 

and derive from Eq. ( ^0|) 

doti , A / 1 i A \ ~ 

-E|f fii^+E^r 1 ^), (22) 

j,k P r=2 71 

where (a) (r = 2, • • •) are power series of a's and can be computed from the r-th loop 
A functions. Next we search for fixed points pi of Eq. (EI) at a = 0. To this end, we have 
to solve 

(-l + |i))^-E|f PiPk = 0, (23) 

3v^ 



7 



and assume that the fixed points have the form 

Pi = for i = 1, • • • , n' ; pi > for i = n' + 1, ■ ■ • , n . (24) 

We then regard <5j with i < n' as small perturbations to the undisturbed system which is 
defined by setting a,i with i < n' equal to zero. As we have seen, it is possible to verify 
at the one-loop level |T7J the existence of the unique power series solution 



a i = Pi + Pi a> r 1 , i = n' + l,---,n (25) 

r=2 

of the reduction equations fl22|) to all orders in the undisturbed system. These are RGI 
relations among couplings and keep formally perturbative renormalizability of the undis- 
turbed system. So in the undisturbed system there is only one independent coupling, the 
primary coupling a. 

The small perturbations caused by nonvanishing <Sj with i < n' enter in such a way 
that the reduced couplings, i.e., <Sj with i > n', become functions not only of a but also 
of &i with i < n'. It turned out that, to investigate such partially reduced systems, it is 
most convenient to work with the partial differential equations 

~ d n ' ~ d „ „ ~ 
{ P tt + J2 Pa ttz- } &i(a,a) = (3i(a,a) , 
da ^ da a 

n A(o) P ~ n P /r,r\ 

Pi(a) = — I n ' P = - ' ( 26 ) 



which are equivalent to the reduction equations (f22"D, where we let a, b run from 1 to n' 
and i,j from n' + 1 to n in order to avoid confusion. We then look for solutions of the 
form 



ft; 



Pi + Ysi-T' 1 ti T \*») , i = n' + n, (27) 



r=2 * 



where fj; r \a a ) are supposed to be power series of a a . This particular type of solution 
can be motivated by requiring that in the limit of vanishing perturbations we obtain the 
undisturbed solutions (p5|) p9| |29| . Again it is possible to obtain the sufficient conditions 



(r) 

for the uniqueness of /, in terms of the lowest order coefficients. 
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3 Finite Unified Model Based on 577(5) 

According to the discussion in the previous chapter, the non-renormalization theorem en- 
sures that there are no extra mass and cubic-interaction-term renormalizations, implying 
that the /3-functions of can be expressed as linear combinations of the anomalous 
dimension matrix 7^ of </>*. Therefore, all the one-loop /3-functions of the theory vanish if 

/?« = and 7 « = (28) 



are satisfied, where f3^ and 7^ are given in Eqs. (0) and ([18) respectively. A very 
interesting result is that these conditions (|28| ) are necessary and sufficient for finiteness 



at the two-loop level p2] 



A natural question is what happens in higher loops. Since the finiteness conditions 
impose relations among couplings, they have to be consistent with the REs (^]) (this should 
be so even for the one-loop finiteness). Interestingly, there exists a powerful theorem |27[ 
which provides the necessary and sufficient conditions for finiteness to all loops. The 
theorem makes heavy use of the non-renormalization property of the supercurrent anomaly 



28]. In fact, the finiteness theorem can be formulated in terms of one-loop quantities, and 



it states that for supersymmetry gauge theories we are considering here, the necessary 



and sufficient conditions for (3 g and (3ijk to vanish to all orders are [27| : 

(a) The validity of the one-loop finiteness conditions, i.e., Eq. (^) is satisfied. 

(b) The REs (2) admit a unique power series solution, i.e., the corresponding matrix M 
defined in Eq. (|9]) with (3^ = has to be non-singular. 

The latter condition is equivalent to the requirement that the one-loop solutions p^'s are 
isolated and non-degenerate. Then each of these solutions can be extended, by a recursion 
formula, to a formal power series in g giving a theory which depends on a single coupling 
g, and has /3-functions vanishing to all orders. 

From the classification of theories with M 1 ) = [23], one can see that using SU(5) as 



gauge group there exist only two candidate models which can accommodate three fermion 
generations. These models contain the chiral supermutiplets 5 , 5 , 10 , 5 , 24 with the 
multiplicities (6, 9, 4, 1, 0) and (4, 7, 3, 0, 1), respectively. Only the second one contains a 
24-plet which can be used for spontaneous symmetry breaking (SSB) of SU(5) down to 
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£77(3) x SU(2) x £7(1). (For the first model one has to incorporate another way, such as 
the Wilson flux breaking to achieve the desired SSB of £77(5).) Here we would like to 
concentrate only on the second model. 

The most general £77 (5) invariant, cubic superpotential of the (second) model is: 

W = H a [f ab H b 24 + h ia 5a4 + g ija 10,5,] +p(24) 3 

+ - 10; [ g ija lOjHa + g iab Ejl h + g' ijk 5j5 k ] , (29) 

where i,j, k = 1, 2, 3 and a, b = 1, • • • , 4, and we sum over all indices in W (the SU (5) 
indices are suppressed). The 107s and 57s are the usual three generations, and the four 
(5 + 5) Higgses are denoted by H a , H a . 

Given the superpotential, the 7^ 1 ' ), s can be easily computed {(3^ vanishes of course). 
To ensure finiteness of the model to all orders, we have to find pW's that are isolated and 
non-degenerate solutions of Eq. (6) and are consistent with the vanishing 7 ( 7)'s. In most 
of the previous studies of the present model [24", 25], however, no attempt was made to 



find isolated and non-degenerate solutions, but rather the opposite. They have used the 
freedom offered by the degeneracy in order to make specific ansatze that could lead to phe- 
nomenologically acceptable predictions. Here we concentrate on finding an isolated and 
non-degenerate solution that is phenomenologically interesting. As a first approximation 
to the Yukawa matrices, a diagonal solution, that is, without intergenerational mixing, 
may be considered. It has turned out that this can be achieved by imposing the Z 7 x Z% 
discrete symmetry and a multiplicative Q-parity on W, and that, in order to respect these 
symmetries, only gm , g iU , fa and p are allowed to be non-vanishing. Moreover, we have 
found that under this situation there exists a unique reduction solution that satisfies the 
finiteness conditions (a) and (b) J7J: 

8 6 
ocui = -atauT + 0(a 2 GUT ) , am = -a GUT + 0(a 2 GUT ) , a hi = , 

15 

a/44 = a GUT + 0(a GUT ) , a p = —a GUT + 0(a 2 GUT ) , (30) 

where i = 1,2,3, and the 0(a GUT ) terms are power series in cxgut that can be uniquely 
computed to any finite order if the /3-functions of the unreduced model are known to 
the corresponding order. The reduced model in which gauge and Yukawa couplings are 
unified has the /3-functions that identically vanish to that order. 
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In the above model, we found a diagonal solution for the Yukawa couplings, with each 
family coupled to a different Higgs. However, we may use the fact that mass terms do 
not influence the /3-functions in a certain class of renormalization schemes, and introduce 
appropriate mass terms that permit us to perform a rotation in the Higgs sector such that 
only one pair of Higgs doublets, coupled to the third family, remains light and acquires a 
non- vanishing VEV (in a similar way to what was done by Leon et al. ||25|| ). Note that the 
effective coupling of the Higgs doublets to the first family after the rotation is very small 



avoiding in this way a potential problem with the proton lifetime Thus, effectively, we 
have at low energies the MSSM with only one pair of Higgs doublets. Adding soft breaking 
terms (which are supposed not to influence the /3-functions beyond M GUT ) , we can obtain 
supersymmetry breaking. The conditions on the soft breaking terms to preserve one-loop 



finiteness have been given already some time ago ||32| . Recently, the same problem at 
the two- loop level has been addressed |33|]. It is an open problem whether there exists 
a suitable set of conditions on the soft terms for all-loop finiteness. Since the £77(5) 
symmetry is spontaneously broken below Mgut, the finiteness conditions obviously do 
not restrict the renormalization property at low energies, and all it remains is a boundary 
condition on the gauge and Yukawa couplings; these couplings at low energies have to be 
so chosen that they satisfy fl30l) at M GUT . So we examine the evolution of the gauge and 
Yukawa couplings according to their renormalization group equations at two-loops taking 
into account all the boundary conditions at Mqut- The predictions for mt for a varying, 
but common to all particles, supersymmetry breaking threshold are given in Figure 1. 



4 The Minimal Asymptotically Free SU(5) Model 

Let us consider next the minimal N = 1 supersymmetric gauge model based on the group 
££7(5) ||. Its particle content is then specified and has the following transformation prop- 
erties under £77(5): three (5 + 10)- supermultiplets which accommodate three fermion 
families, one (5 + 5) to describe the two Higgs supermultiplets appropriate for electroweak 
symmetry breaking and a 24-supermultiplet required to provide the spontaneous symme- 
try breaking of SU(5) down to SU(3) x SU(2) x 17(1). 
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Since we are neglecting the dimensional parameters and the Yukawa couplings of the 
first two generations, the superpotential of the model is exactly given by 

W = X -g t 10 3 10 3 H + g b 5 3 10 3 H + g x (24) 3 + g f H 24 H , (31) 

where H, H are the 5, 5- Higgs supermultiplets and we have suppressed the SU(5) indices. 

According to the notation introduced in Eq. (JZT]), Eqs. (^) become 

da t 27 _ _ 2 4 _ _ 8 _ _ 
a — — = — a* — 3 a t — - a t a b — - a t at , 
da 5 3 5 

da 6 23 _ 10 _ 2 _ 8 „ _ 

a — — = — ah a h — a b at ah at , 

da 5 3 b 5 f ' 

n ~ 21 ^ 2 _ _ 
a—— = 9a x - — a x - a x a f , 
da 5 

datf 83 „ 53 ~ 7__ . . 

""cfcT = 15 a/ _ 15 af ~ af0tt ~ 3 af0Lb ~ 5 a/ " A ' 

in the one-loop approximation. Given the above equations describing the evolution of 

the four independent couplings («j , i — t,b,X,f), there exist 2 4 = 16 non-degenerate 

solutions corresponding to vanishing p's as well as non-vanishing ones given by Eq. (|27|) 

The possibility to predict the top quark mass depends, as in the previous model, on 

a nontrivial interplay between the vacuum expectation value of the two SU(2) Higgs 

doublets involved in the model and the known masses of the third generation (m& , m T ). 

It is clear that only the solutions of the form 

Pt , Pb ± (33) 

can predict the top and bottom quark masses. 

There exist exactly four such solutions. The first solution is ruled out since it is 
inconsistent with Eq. ([19]), and the second one is ruled out since it does not satisfy the 
criteria to be asymptotically free. We are left with two asymptotically free solutions, 
which we label 3 and 4. According to the criteria of section 2, these two solutions give 
the possibility to obtain partial reductions. To achieve this, we look for solutions || of 
the form Eq. (|25|) to both 3 and 4. 

We present now the computation of some lower order terms within the one-loop ap- 
proximation for the solutions. For solution 3: 

5i = Vi + ft x=1) &x + ft x=2) &l + --- foii = t,b,f, (34) 
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where 



Vt,b,f 



2533 1491 560 



2605 ' 2605 ' 521 ' 



fi r b x J l) ^ 0.018,0.012, -0.131, 

4?/ =2) - °- 005 > °- 004 > - °- 021 > (35) 



For the solution 4, 



where 



&i = Vi + ft f=1) a f + ft> =1) ax + fp =1 ^ =1) a, a x 

+ft f=2) <% + ft x=2) *i--- far < = *,&, (36) 



89 63 (rA= i) {rx=2) 

'lt,b nr i nr j Ji Ji u i 

65 65 

// r b /=1) ~ -0.258 , - 0.213 , fll 1=l) ~ - 0.258 , - 0.213 , 

//; /=2) ~ -0.055 , - 0.050 , ft r b f=1 ' rx=1) ^ - 0.021 , - 0.018 , (37) 



In the solutions ( p4[) and (^g) we have suppressed the contributions from the Yukawa 
couplings of the first two generations because they are negligibly small. 

Presumably, both solutions are related; a numerical analysis on the solutions || sug- 
gests that the solution 3 is a "boundary" of 4. If it is really so, then there is only one 
unique reduction solution in the minimal supersymmetric GUT that provides us with the 
possibility of predicting a t . Note furthermore that not only a t but also af, is predicted in 
this reduction solution. 

Just below the unification scale we would like to obtain the standard SU(3) x SU(2) x 
U(l) model while assuming that all the superpartners are decoupled at the Fermi scale. 
Then the standard model should be spontaneously broken down to SU(3) x U(l) em due 
to VEV of the two Higgs S , f/(2)-doublets contained in the 5, 5-super-multiplets. One way 
to obtain the correct low energy theory is to add to the Lagrangian soft supersymmetry 
breaking terms and to arrange the mass parameters in the superpotential along with the 
soft breaking terms so that the desired symmetry breaking pattern of the original SU(5) is 
really the preferred one, all the superpartners are unobservable at present energies, there 
is no contradiction with proton decay, and so forth. Then, as in the previous model, we 
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study the evolution of the couplings at two loops respecting all the boundary conditions 
at Mqut- The predictions for m t versus Msusy for the two sets of boundary conditions 
given above (AFUT3 and AFUT4) together with the corresponding predictions of the 
FUT model, are given in Figure 1. In a recent study ||34|| , we have considered the proton 



decay constraint [35[] to further reduce the parameter space of the model. It has been 
found that the model consistent with the non-observation of the proton decay should be 
very close to AFUT3, implying a better possibility to discriminate between the FUT and 
AFUT models, as one can see from Figure 4. 



5 Asymptotically Non-Free Supersymmetric Pati- 
Salam Model 

In order for the RGI method for the gauge coupling unification to work, the gauge cou- 
plings should have the same asymptotic behavior. Note that this common behavior is 
absent in the standard model with three families. A way to achieve a common asymp- 
totic behavior of all the different gauge couplings is to embed SU(3)c x SU{2) L x U{l)y 
to some non-abelian gauge group, as it was done in sections 3 and 4. However, in this 
case still a major role in the GYU is due to the group theoretical aspects of the covering 
GUT. Here we would like to examine the power of RGI method by considering theories 
without covering GUTs. We found |9| that the minimal phenomenologically viable model 
is based on the gauge group of Pati and Salam Q PS = SU(4) x SU{2) R x SU{2) L . We 
recall that N = 1 supersymmetric models based on this gauge group have been studied 
with renewed interest because they could in principle be derived from superstring [3~7|. 



In our supersymmetric, Gauge- Yukawa unified model based on Q PS 0, three gener- 
ations of quarks and leptons are accommodated by six chiral supermultiplets, three in 
(4, 2, 1) and three (4, 1, 2), which we denote by vJ/Wm lR and ^>^ %L . (I runs over the three 
generations, and /i, v (= 1,2,3,4) are the S77(4) indices while i R , ii (= 1,2) stand for 
the SU(2) LjR indices.) The Higgs supermultiplets in (4, 2, 1), (4, 2, 1) and (15, 1, 1) are 
denoted by lR , iR and S^, respectively. They are responsible for the spontaneous 
symmetry breaking (SSB) of SU(4) x SU(2) R down to SU(3) C x U{l) Y . The SSB of 
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U(1)y x SU(2)l is then achieved by the nonzero VEV of hi R i L which is in (1,2,2). In 
addition to these Higgs supermultiplets, we introduce G^ i iL (15, 2, 2) , <fi (1, 1, 1) and 
E'jf (15, 1, 1). The Gl iRih is introduced to realize the SU{4) x SU{2) R x SU{2) L version 



of the Georgi-Jarlskog type ansatz [38] for the mass matrix of leptons and quarks while 
is supposed to mix with the right-handed neutrino supermultiplets at a high energy scale. 
With these in mind, we write down the superpotential of the model W, which is the sum 
of the following superpotentials: 

3 



W Y = £ gjJ ¥»* h iRiL , W GJ = g GJ ^f lR G$ iRjL *W» * 

i,j=i 

(I)iR 



W NM = £ gi* s * 

7=1,2,3 

— Tr [ E 3 1 + — ' r|V 1 ' v '' r - 
3 L J 2 



W SB = gHH^^H^v + ^Tr [ S 3 ] + ^ Tr [ (E') 2 E 



W TDS = ^ e iR3R e iL3L Xr [ G iRiL E G jRjL ] , 
V^m = m h h 2 + m G G 2 + m^<f) 2 + m H HH + mzZ 2 + m J ; l {i:') 2 . (38) 

Although has the parity, <fi — > —<ft and E' — > — S', it is not the most general potential, 
and, by virtue of the non-renormalization theorem, this does not contradict the philosophy 
of the coupling unification by the RGI method. 

We denote the gauge couplings of SU(4) x SU(2)r x SU{2)l by a.\ , olir and 0:21, 
respectively. The gauge coupling for Z7(l)y, a%, normalized in the usual GUT inspired 
manner, is given by 1/aii = 2/5a4 + 3/5a2R ■ In principle, the primary coupling can 
be any one of the couplings. But it is more convenient to choose a gauge coupling as 
the primary one because the one-loop /3 functions for a gauge coupling depends only on 
its own gauge coupling. For the present model, we use «2L as the primary one. Since 
the gauge sector for the one-loop j3 functions is closed, the solutions of the fixed point 
equations fl2"3D are independent on the Yukawa and Higgs couplings. One easily obtains 
Pa* 1 = 8/9 , P2R — 4/5, so that the RGI relations ( p7|) at the one- loop level become 

"48. ai 5 
"4 = = r , "i = = -x ■ (39) 

The solutions in the Yukawa-Higgs sector strongly depend on the result of the gauge 
sector. After slightly involved algebraic computations, one finds that most predictive 
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solutions contain at least three vanishing pf^ 's. Out of these solutions, there are two that 
exhibit the most predictive power and moreover they satisfy the neutrino mass relation 



m,. 



> m L 



m. 



For the first solution we have p\ 



(i) 



0, while for the 



second solution, p^) = p^J, = Pq j = 0, and one finds that for the cases above the power 
series solutions take the form 



«33 



"30 



1.67 - 0.05a^ + O.OO45 2 - 0.905 E + ■ 
2.20 - O.O8a 2 - 0.055 G H 

3.33 + O.O5«i + O.21a 2 - 0.025 s + • • ■ 
3.40 + O.O55 10 - 1.63520 - 0.0015 G + • 

1.43 - 0.585^ - 1.435 2 - 0.035 2 + ■ ■ 
0.88 - 0.485^ + 8.835 20 + 0.015 G + ■ • 

1.08 - 0.035^ + O.1O5 20 - 0.075 s + • ■ • 
2.51 - 0.045^ - 1.68520 - 0.125 G + ■ ■ 

0.40 + 0.015^ - O.455 20 - 0.105 G + • • 

4.91 - 0.0015^ - O.O35 2 - 0.465 s + ■ 
8.30 + 0.015x0 + 1.725 2 - 0.365 G + • • 

5.59 + O.O25i0 - O.O45 2 - 1.335 s + • • 



(40) 



We have assumed that the Yukawa couplings gjj except for g 33 vanish. They can be 
included into RGI relations as small perturbations, but their numerical effects will be 
rather small. 

The number Njj of the Higgses lighter than Msusy could vary from one to four while 
the number of those to be taken into account above Msusy is fixed at four. We have 
assumed here that Nh = 1. The dependence of the top mass on Msusy in this model is 
shown in Figure 2. 
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6 Discussion and Conclusions 



As a natural extension of the unification of gauge couplings provided by all GUTs and 
the model dependent unification of Yukawa couplings, we have introduced in a number 
of publications the idea of Gauge- Yukawa Unification. GYU is a functional relationship 
among the gauge and Yukawa couplings provided by some principle. In our studies GYU 
has been achieved by demanding the principles of reduction of couplings and finiteness. 
The first principle is based on the existence of RGI relations among couplings which 
preserve perturbative renormalizability in gauge theories. The second principle is based 
on the fact that it is possible to find RGI relations among couplings that keep finiteness 
in perturbation theory, even to all orders. In the previous chapters we have presented the 
application of these principles in various models as well as the resulting predictions. 

The consequence of GYU is that in the lowest order in perturbation theory the gauge 
and Yukawa couplings above Mqut are related in the form 

9i = lugaur , i = l,2,3,e, • • -,r,b,t , (41) 

where gi (i = 1, •••,£) stand for the gauge and Yukawa couplings, gcuT is the unified 
coupling, and we have neglected the Cabibbo-Kobayashi-Maskawa mixing of the quarks. 
So, Eq. (f41~|) exhibits a boundary condition on the the renormalization group evolution 
for the effective theory below Mqut; which we assume to be the MSSM. 

As we have seen in the previous chapters, there are various supersymmetric GUTs 
with GYU in the third generation that can predict the bottom and top quark masses 
in accordance with the experimental data [J. This means that the top-bottom 

hierarchy could be explained in these models, exactly in the same way as the hierarchy of 
the gauge couplings of the SM can be explained if one assumes the existence of a unifying 
gauge symmetry at Mqut @- 

It has been also observed [|l(J that there exists a relatively wide range of k's which 
gives the top-bottom hierarchy of the right order. Of course, the existence of this range 



is partially related to the infrared behavior of the Yukawa couplings [35]. Therefore, a 
systematic investigation on the nature of GYU is indispensable to see whether a GYU 
can make experimentally distinguishable predictions on the top and bottom masses, or 
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whether the top-bottom hierarchy results mainly from the infrared behavior of the Yukawa 
couplings. With more precise measurements of the top and bottom masses, we will be 
able to conclude which case is indeed realized. 



We have performed an exhaustive analysis on this problem at the two-loop level [10 



and here we would like to present only a few representative results to provide an idea 
of our complete analysis. We have made the same assumptions as in chapters 3, 4, and 
5, namely that below Mqut the evolution of couplings is governed by the MSSM and 
that there exists a unique threshold M SU sy for all superpartners of the MSSM so that 
below MgusY the SM is the correct effective theory, where we include only the logarithmic 
and two-loop corrections for the RG evolution of couplings. We have neglected all the 
threshold effects. Note that with a GYU boundary condition alone the value of tan (3 can 
not be determined; usually, it is determined in the Higgs sector, which however strongly 
depends on the supersymmetry breaking terms. In our analysis we avoid this by using 
the tau mass, along with M z , a~^(M z ) and sin 2 9 w (Mz), as the input. Figure 3 shows 
the dependence of the top mass on different values of k 2 , when n 2 = k 2 /k 2 * s fixed to 2.0, 
and the supersymmetry breaking scale Msusy = 500 GeV. 

At this point it is also worth recalling the predictions for m t of ordinary GUTs, in 
particular of supersymmetric SU(5) and 50(10). The MSSM with SU(5) Yukawa bound- 
ary unification allows m t to be anywhere in the interval between 100-200 GeV for 
varying tan/5, which is now a free parameter. Similarly, the MSSM with 5*0(10) Yukawa 
boundary conditions, i.e. t — b — r Yukawa Unification gives m t in the interval 160-200 
GeV ||. 

It is clear that the GYU scenario presented here, provides us with predictive solutions 
which are consistent with all relevant experimental data. Moreover, it is the most predic- 
tive scheme as far as the mass of the top quark is concerned. It is worth noting that the 
various GYU models can be well distinquished among themselves when more information 
on the supersymmetry breaking scale is available. A nice demonstration of this point is 
shown in Figure 4. 

Concerning recent related studies, we would like to emphasize that our approach of 
dealing with asymptotically non-free theories || covers work done by other authors f4lf , 
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though the underlying idea might be different. Also, interesting RGI relations among the 



soft breaking parameters above the unification scale have been found |42|]. These relations 
are obtained in close analogy to our approach. 

It will be very interesting to find out in the comming years, as the experimental 
accuracy of m t increases, if nature is kind enough to verify our conjectured Gauge- Yukawa 
Unification. 
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Figure Captions 



Fig. 1. The comparison of the predictions for m t between the AFUT models and the 
FUT one. For the FUT model a t = 1.6, a b = 1.2, for AFUT3 a t = 0.97, a b = 0.57, and 
for AFUT4 a t = 1.37, & b = 0.97. 

Fig. 2. The values for m t predicted by the Pati-Salam model for different M S usy scales. 
Only the ones with M SUSY beyond 400 GeV are realistic. 

Fig. 3. The dependence of the top mass m t with k\ , at fixed M S usy = 500 GeV. As we 
can see, after k\ ~ 2.0 the top mass goes to its infrared fixed point value. 

Fig. 4. m t predictions of SU(5) FUT and AFUT3 models, for given M SUS y around 
100 and 500 GeV. For the FUT model a t = 1.6, a b = 1.2, and for AFUT3 a t = 0.97, 
a b = 0.57. 
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